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Abstract 

We have studied the structure and properties of potassium clusters containing even number 



of atoms ranging from 2 to 20 at the ab initio level. The geometry optimization calculations are 
q performed using all-electron density functional theory with gradient corrected exchange-correlation 

•i-H 1 

^afunctional. Using these optimized geometries we investigate the evolution of binding energy, ioniza- 
Qjtion potential, and static polarizability with the increasing size of the clusters. The polarizabilities 
are calculated by employing Moller-Plesset perturbation theory and time dependent density func- 
J_. tional theory. The polarizabilities of dimer and tetramer are also calculated by employing large 
£vq basis set coupled cluster theory with single and double excitations and perturbative triple exci- 
tations. The time dependent density functional theory calculations of polarizabilities are carried 
©ut with two different exchange-correlation potentials: (i) an asymptotically correct model po- 

tential and (ii) within the local density approximation. A systematic comparison with the other 

> ■ 

Available theoretical and experimental data for various properties of small potassium clusters men- 
C$ tioned above has been performed. These comparisons reveal that both the binding energy and the 
ionization potential obtained with gradient corrected potential match quite well with the already 
published data. Similarly, the polarizabilities obtained with Moller-Plesset perturbation theory 
and with model potential are quite close to each other and also close to experimental data. 
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I. INTRODUCTION 



During last two decades rapid progress in the experimental methods of producing atomic 
and molecular clusters in controlled fashion along with development of sophisticated theo- 
retical tools to handle such finite fermionic systems at the ab initio level led to emergence of 
the field of cluster science as one of the most exciting and productive discipline of physics, 
chemistry, and material science. Metal clusters, specially those of alkali-metal atoms Li, 
Na, and K played an important role in the development of cluster physics as a branch of 
modern physics and chemistry. Interest in the study of alkali-metal clusters grew with the 
pioneering work of Knight and co-workers |l| . These researchers discovered that certain 
clusters, those with magic number 8, 20, 34, 40, ■ • ■ of atoms are more stable and conse- 
quently were found more abundantly in the mass spectra of these clusters. The existence of 
magic number clusters is attributed to electronic shell structure of the clusters. Other prop- 
erties like ionization potential, electron affinity, and static polarizabilities of metal clusters 
also show significance of the shell structure. Besides these, photoabsorption cross sections 
have also been measured for alkali metal clusters and have been investigated theoretically 
at various levels. A large body of theoretical work on the electronic structure and optical 
response properties of alkali-metal clusters exists in the literature. Majority of the theoret- 
ical work have been carried out by employing Density functional theory (DFT) and time 
dependent DFT (TDDFT) within the spherical jellium background model (SJBM) (see the 
review article j^, 3]). The SJBM replaces the discrete ionic structure of clusters by a spher- 
ically symmetric uniform positive charge background thus making it possible to carry out 
calculations for the optical response properties of reasonably large clusters of around 100 
atoms 3, 4]. In last ten years or so, several all electron ab initio calculations devoted to the 



ground state and the optical response properties of sodium clusters taking into account the 
actual geometrical arrangement of the sodium atoms have been reported in the literature 

b, a i, a y y y & y a ^ ^ cou , d hmd , ^ ^ 

sizes smaller than those that could be studied by performing jellium based calculations. 
We should mention here that the calculations of structure, electronic and static polraizabil- 
ity of small sodium and lithium clusters have also been performed by employing; ab initio 



correlated wavefunction based methods like MP2, MP4, and SCF-CI 
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We note here that among these numerous studies involving properties of alkali-metal 
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clusters only very few are devoted to the calculations of properties of potassium atom clus- 
ters. This is quite surprising considering the fact that the experimental results for the 
ionization potential, static polarizability, and photoabsorption spectra of potassium clus- 
ters as functions of cluster size were reported very early in the development of cluster 
physics Along with the calculations of sodium clusters few SJBM based stud- 

ies within the realm of DFT and TDDFT pertaining to the evolution of binding ener- 
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s, static polarizabilities of potassium clusters exist in the literature 
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261 ] . Likewise, few studies of the potassium clusters at the ab initio 



level taking into account the detailed ionic structure employing various approaches like con- 



figuration interaction (CI) 27], |28j], many-body perturbation theory 29], DFT coupled with 



pseudopotential [30] , and coupled cluster theory [3JJ have also been reported in the literature 
. However, these studies were restricted to very small sized clusters containing maximum up 
to seven potassium atoms [29J] and investigated ground state properties like bond lengths, 
binding energies, and ionization potentials. To the best of our knowledge no ab initio level 
calculation of static polarizability either by employing correlated wavefunction approach or 
by DFT/TDDFT methods exists in the literature. The aim of this paper is to fill this gap 
by extending the ab initio study to larger potassium clusters containing up to 20 potassium 
atoms and investigate systematically the size evolution of various ground state properties 
like the binding energy (BE), ionization potential (IP) and response property like static po- 
larizability. For this purpose we employ ab initio DFT based method with gradient corrected 
exchange-correlation (XC) potential for geometry optimization and three methods for the 
calculations of polarizabilities namely, the second order Moller-Plesset perturbation theory 
(MP2) [32|] , coupled cluster theory with single and double excitations and perturbative triple 
excitations (CCSD(T)) 33J, and TDDFT with asymptotically correct XC potential. The 



high computational cost and time involved in CCSD(T) and MP2 based calculations restrict 
the maximum size of the cluster that could be handled by this approach. In this paper we 
perform MP2 based polarizability calculations for clusters containing maximum up to 14 
atoms and CCSD(T) based calculations could be performed for dimer and tetramer only. 
Nonetheless these results provide a way to check the correctness and consistencies of the 
results obtained by us by employing various methods. 

We note that to carry out ab initio calculations of ground state properties mentioned 
above and the polarizabilities of clusters which go beyond jellium model, it is necessary 
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to have the knowledge of the ionic structure of the clusters. To this end, we carry out 
a systematic search for the optimized structures by employing DFT based geometry opti- 
mization scheme with Becke-Perdew (BP86) 34| exchange-correlation (XC) potential within 
the generalized gradient approximation (GGA) as this potential is known to yield reliable 
geometries. For each cluster beyond K 4 several structures have been considered as start- 
ing geometries. To the best of our knowledge studies involving geometry optimization of 
potassium clusters beyond K7 are not available in the literature and our results in this pa- 
per will be useful for further studies on potassium clusters. A brief discussion on geometry 
optimization procedure is given in next section. 

The calculation of response property like polarizabilities by TDDFT approach requires 
approximating the forms of XC functionals. It is well known that the accuracy of the results 
for response properties obtained via TDDFT crucially depends on the nature of the XC 



potential, especially its behaviour in the asymptotic region [35|, |36( . Keeping this in mind, 
we carry out all-electron TDDFT based calculations of static poarizabilities of potassium 
clusters with a model XC potential, called statistical average of orbital potentials (SAOP), 
which has desirable properties both in the asymptotic and the inner regions of a molecule 



23, 



38] . The choice of SAOP is also motivated by the results of Ref . lOj, ll5| where it has been 



shown that SAOP polarizability and excitation energies of sodium clusters agree well with 
the experimental data. In order to study the effect of XC potential on the polarizabilities, 
calculations of polarizabilities are also carried out with less accurate XC potential under 
local density approximation (LDA). 

It is well known that the static polarizabilities of sodium clusters obtained by employing 
DFT and TDDFT within SJBM are generally underestimated in comparison to the corre- 
sponding ah initio and experimental results. However, such comparison of the SJBM, ah 
initio, and experimental results does not exist for potassium clusters. In order to test the 
accuracy of the jellium model, we compare the results for the static polarizability of 8 and 
20 atom clusters (as jellium based results only for these two magic clusters are available 
in the literature) obtained by employing LDA XC potential in the realm of SJBM with 
corresponding ah initio and experimental data. 

The rest of the paper is organized as follows: In section II we discuss the theoretical 
methods employed to calculate the optimized geometry, BE, IP, and static polarizabilities 
of potassium clusters. Results of our calculations are presented in Section III and the paper 
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is concluded in Section IV. 



II. THEORETICAL METHODS 

In this work we study the structure and other properties of small sized potassium clusters 
on the basis of all-electron ab initio methods. For this purpose we use ADF program package 



391 ] for DFT and TDDFT based calculations and GAMESS electronic structure code for 



carrying out post-Hartree-Fock MP2 and CCSD(T) calculations 40|. We calculate the 
optimized geometries of clusters containing even number of atoms ranging from 2 to 20 
atoms. The geometry optimizations of all the clusters have been performed through DFT 
based calculations by employing a triple-^ Slater-type orbital (STO) basis set with two added 
polarization functions (TZ2P basis set of ADF basis set library) along with Becke-Perdew 
(BP86) XC potential [34|. The geometry optimization involves finding local minima on the 
multidimensional potential energy surface. The staring geometry of the cluster plays an 
important role in the optimization procedure. In the present calculations for each cluster 
we have considered more than one starting geometries which are compiled from already 



available optimized structures of sodium clusters [11|, [12j. Here we make no assumption 
regarding the core electrons and perform all-electron geometry optimization calculations. 
All the optimizations are carried out with the convergence criteria for the norm of energy 
gradient and energy, fixed at 10~ 4 atomic units (a.u.) and 10~ 6 a.u., respectively. 

In this paper we further calculate the static polarizabilities of potassium clusters by em- 
ploying MP2, CCSD(T), and TDDFT based methods. The calculations of polarizabilities 
with MP2 and CCSD(T) have been carried out by employing finite field approach available 
in the GAMESS electronic structure code. The finite field approach makes use of the per- 
turbative series expansion of the energy E in terms of the components of a static uniform 
electric field F given by, 

E{F) = E(0) + ^ + \Y1 a ^ F 3 + ■ ■ • > (!) 

i ij 

where E(0) is the energy of the system in the absence of the applied electric field, jl is dipole 
moment, and a^- (i,j = x,y,z) is the dipole polarizability tensor. The components of the 
polarizabilty tensor are obtained as the second-order derivatives of the energy with respect 
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to the components of the electric field, 

The derivatives are calculated numerically by applying fields of 0, 0.001 and 0.002 a.u. along 
±x, ±y, and ±z directions and mean polarizability is calculated from the diagonal elements 
of polarizability tensor as 



« = 77 + oi yy + a zz ) . (3) 



1 

3 

All MP2 and CCSD(T) calculations have been performed with valence triple zeta polarized 



Gaussian basis sets of Sadlej and Urban [41] for potassium atom. 

On the other hand, TDDFT calculation of polarizability is based on the linear response 
theory of many-body systems and employs exact analytical expressions for polarizability 
in terms of the moment of the first-order induced density. To avoid digression we refer the 
readers to Ref. [42| for detailed description of the linear response theory based method which 
is adopted in ADF program package for obtaining polarizability. We note here that TDDFT 
based method gives frequency dependent polarizability but here we focus our attention on 
the static or zero-frequency polarizability. We also mention here that the calculation of 
polarizability based on linear response theory of many-body system yields results which are 
more accurate than the finite field approach as no explicit specification of the magnitude 
of the applied field is required in the former method. It has been already mentioned that 
TDDFT based response property calculation requires approximating the XC functional at 
two different levels. The first one is the static XC potential needed to calculate the ground- 
state Kohn-Sham (KS) orbitals and orbital energies. The second approximation is needed to 
represent the XC kernel fxc{ v i r u ) which determines the XC contribution to the screening 
of an applied field. For the XC kernel, we use reasonably accurate adiabatic local density 



approximation (ALDA) 43j. On the other hand, for the static XC potential needed to 



calculate the ground-state orbitals and energies, two different choices have been made. These 
are (i) the standard potential under local density approximation (LDA) as parametrized 
by Vosko, Wilk and Nussair [44] and (ii) the model potential SAOP possessing correct 



23, 



381 ] . The results obtained by these two 



behaviour both inner and asymptotic regions 
XC potentials are compared in order to investigate the effect of XC potential on the results 
for the polarizability. The calculations of polarizabilities of potassium clusters by TDDFT 
based method are carried out by using large Slater type orbital (STO) basis sets. It is well 



known that for accurate calculations of response properties it is necessary to have large 
basis sets with both polarization and diffuse functions. For our purpose, we have chosen one 
of the largest all electron even tempered basis set ET-QZ3P-2DIFFUSE with two sets of 
diffuse functions consisting of (13s,10p,5d,3f) functions for K available in the ADF basis set 
library. The application of basis set with diffuse functions often leads to the problem of linear 
dependencies. Such problems have been circumvented by removing linear combinations of 
functions corresponding to small eigenvalues of the overlap matrix. We expect that the 
size of the chosen basis set will make our results very close to the basis-set limit. The 
next section is devoted to the discussion of the results for the structures and properties of 
potassium clusters. 



III. RESULTS AND DISCUSSION 



We begin this section with the discussion on the results of our geometry optimization 
calculations followed by the results for BE, IP, and polarizabilities of clusters consisting of 
even number of atoms ranging from 2 to 20. We compare the results of our calculations with 
the available experimental data and also results of other theoretical works performed both 
within the framework of the jellium model and beyond, using DFT or correlated wavefunc- 
tion based methods, and assess the level of accuracy of different theoretical approaches. 



A. Structure and other properties of K n clusters 

The optimized structures of potassium clusters with even number of atoms up to 20 are 
shown in Fig. 1 (n =2, 4, and 6) and Fig. 2 (n > 8). The indices n and m in the label 
n_m assigned to each cluster in these figures denote the number of atoms and the rank in 
the increasing energy order. In order to perform the geometry optimization calculations of 
potassium clusters, we make use of the optimized geometries of sodium clusters reported in 
Refs. I 111 . 112 1 , as initial geometries. For each cluster, we consider all the isomers available 



in Refs. 



nl, 



12j | to search for the minimum energy structures. We note here that structures 
reported in these two papers for sodium clusters are quite exhaustive and include structures 
obtained by several other authors also. Therefore we feel that the minimum energy structures 
obtained by us represent true ground state structures of the potassium clusters. Now we 
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present the results for BE per atom and the average inter-atomic distances for all the isomers 



for each cluster. We tabulate the results in two parts. In Table I, we present t 



27 



28 



re results for 
30j and 



29 



K 2 , K 4 , and K 6 clusters and compare them with the other theoretical 
experimental (for dimer) 4|j] data available in literature. On the other hand, the results 
for the K 8 and higher clusters are presented in Table II for which no data are avaiable for 
comparison. We calculate BE which is presented in Table I and II by employing the formula 



E h {K n ) = nE(K) - E(K n ) (4) 

where E(K n ) and E(K) are the total energies of a neutral n-atom potassium cluster K n 
and an isolated single potassium atom, respectively. Note that according to EqJH BE is a 
negative number for a bound structure and larger value of it implies a more stable structure. 
The second quantity for which results are presented in Table I and II is the average inter- 
atomic distance in each isomer of a cluster. This quantity is computed by employing the 
corresponding optimized structure and considering only inter-atomic distances smaller than 
5.076A which is 10% higher than the nearest neighbour disatance in the bcc lattice of bulk 
potassium. 

First we discuss the results for small clusters K 2 , K4, and K$ (Table I) as for these systems 



we can assess the accuracy of our DFT based results by comparing them with t 
data that already exist in the literature for these three clusters. In Refs. 127 



re published 



28, \m 



the calculations on above mentioned clusters were carried out by emplo y in g pseudopotential 
method in conjunction with configuration interaction (CI) approach |27l. |28| . self interaction 
corrected DFT [30j method, and CCSD(T) for the valence electrons [3l|. On the other hand, 
an all-electron calculation using the techniques of Hartree-Fock theory followed by many- 



body perturbation theory (MBP' 
of potassium clusters up to K 7 



) was carried out to determine the equilibrium geometries 



29] . It can be seen from Table I that for dimer, results 



for both BE and average interatomic distance obtained by us employing DFT with BP86 
XC potential are quite close to the other published data including the experimental result 



451 ]. In fact our result for the BE per atom for K 2 matches exactly with the experimental 



data and our result for interatomic distance is closest to the experimental value compared 
to other theoretical data. 



Following previous studies [27J, [28|, [30j , we consider rhomboidal geometry for tetramer K^. 



The average interatomic distance for the tetramer geometry obtained by us is quite close to 



other results and same is true for the BE per atom except for the result of Ref. 27J]. For Kq, 
two geometries, namely a planar and a three-dimensional structures have been considered 
for the geometry optimization. We find that the planar isomer is higher in energy than the 



three-dimensional one in conformity with the results of Refs. |28l . I29j ]. For the minimum 
energy structure of K & the difference between the results for BE and average interatomic 
distance obtained by us and those of Refs. Q, 29] is less than 10%. As mentioned before 
for clusters beyond K7 no results for the optimized geometries of potassium clusters are 
available in the literature for guidance. Consequently, to obtain the geometries of potassium 
clusters containing 8 and more atoms we make use of corresponding geometries of sodium 



clusters reported in Refs. ll|, [12J as the starting geometries and optimize them by employing 



DFT based method with TZ2P basis set and BP86 XC potential. For K 8 we consider two 
stuructures one possessing V>2d symmetry 12] and another having symmetry [ll|. We 



find that D 2 d symmetry isomer of K 8 is lower in energy than structure by around 0.08 
eV. 

For clusters beyond Nas, many (more than two) quasi-degenerate isomers exist. We use 
all these isomers for the geometry optimization of potassium clusters. In order to perform 
the geometry optimization calculations for the clusters beyond K 8 , we sample four isomers 
for K10, three isomers each for K12, K44, Ki6 and Kig and five isomers for K2o- These 
optimized structures in increasing energy order are shown in Fig. 2 and their BE and 
average interatomic distances are presented in Table II. We note here that for K20 the 
higher symmetry T^ isomer does not yield minimum energy and it is for a lower symmetry 
T>2d isomer of 20 atom cluster we get minimum energy. The energy of T^ isomer is around 
0.37 eV higher than the minimum energy D 2f j structure. It is interesting to note that for 
the two magic number clusters Kg and K20, the minimum energy geometries possess D 2 d 
symmetry. 

Having determined the geometry of potassium clusters containing even number of atoms 
up to K20, we next focus our attention on the evolution of ionization potential (IP) with 
the size of clusters, as this dependence has been investigated extensively and these results 



are available in the literature 



46 



47 



481 ] . Apart from these experimental results some 



papers also reported theoretical results for IP of small potassium clusters at the ab initio 
level j^. [30! and also within SBJM 21]. We compare these results with the ones obtained 
in this paper by DFT based calculations with TZ2P basis set and BP86 XC potential. To 
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calculate IP, we restrict ourselves to only the minimum energy isomers for each cluster and 
employ following formula of IP of a cluster containing N atoms: 



IP = E(K n ) + — E(K n ), (5) 

where E(K n ) + and E(K n ) are energies of the singly charged and neutral clusters respectively 
The results of these calculations along with the data available in the literature are presented 
in Table III and Fig. 3. In Table III, we present the results for small sized potassium 
clusters containing 2, 4, and 6 atoms for which both theoretical and experimental results 
are available. For clusters containing more than 6 atoms theoretical results at ab initio level 
are not available, so we compare our results with the experimental data only and these are 
displayed in Fig. 3. From Table III, it can be seen that for K2, K4, and Kg, our DFT based 
results for IP are quite close to other theoretical as well as experimental data. Fig. 3 shows 
;he dependence of cluster IP obtained by us on n along with the experimental results of Ref. 



18 



46 



47 



48( | . This comparison clearly shows that our results follow a similar trend as 
compared to the experimental data with increase in n. However, the theoretical results are 
slightly overestimated with respect to the experimental ones. This may be attributed to the 
fact that the experiments have been performed at finite temperatures [251 ] . Here we wish 
to point out that Fig. 3 also clearly elucidates that IP of potassium clusters decreases with 
increasing cluster size, which is in conformity with the conducting sphere model (CSM) of 



metal cluster 



49 



50 



5l|. To verify this we calculate IP in accordance with the expression 



IP = W + ljl ( fi ) 

where W is the bulk work function and R is radius of the metallic droplet. Following Ref. 
521 ]. we use W = 2.28 eV and R = r^n 1 / 3 ( where = r s is the Wigner-Seitz radius) with 



r s = 4.86a.w.. These results are shown in Fig. 3 by dashed line and our DFT result follow 
a similar trend. 



B. Polarizability of potassium clusters 

In this section we present and discuss the results of our calculations for the static dipole 
polarizabilities of potassium clusters. The static polarizability plays an important role in 
the charctarization of the clusters and it is one of the property which has been extensively 
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54 



551 ] . To the best of our 



measured for metal clusters specially sodium clusters Q, S . 
knowledge the experimental results for static polarizabilities of some potassium clusters are 
available only in Ref. Here we will compare results of our calculations for a of clusters 
containing 2, 8, and 20 atoms with the above-mentioned experimental results. 

We begin our discussion on the results for polarizabilities by first testing the accuracies 
of SAOP and MP2 results for dimer and tetramer against the corresponding data obtained 
with a large basis set coupled cluster theory with single and double excitations and per- 
turbative triple excitation (CCSD(T)) calculations |33j. In Table IV we present the results 
for polarizabilities of K2 and K4 obtained with three different methods along with the other 



theoretical 



31 



56] and experimental 



57l | results for the dimer which are already available 



in literature. From Table IV, we notice that the results for the dimer obtained by CCSD(T) 
approach vary from 486 to 510 a.u. and these results are well within the experimental errors. 
The variation in the different CCSD(T) results is attributed to the use of different basis sets 
and also to the level of calculations (all-electron or pesudopotential). The basis set used 
in the present paper is similar to the one employed in Ref. 56] and both are all-electron 
calculations and consequently two results are close to each other. The difference in the two 
results may be due to the use of slightly different bond length (3.91A) of the dimer and 
also inclusion of relativistic effect in the calculation of Ref. [56]. Furthermore, we observe 
that MP2 value for the polarizability of K 2 is quite close to the lowest CCSD(T) result of 
Ref. 3l|, however, SAOP underestimates the polarizability by around 6% with respect to 
CCSD(T) rseults. In contrast to the dimer case SAOP polarizability for K 4 is quite close to 
the CCSD(T) result and slightly higher (around 2%) than MP2 result. 

We now proceed with the calculations of static polarizabilities for all the optimized iso- 
mers shown in Figs. 1 and 2. We note here that the geometries considered in this paper 
are nonspherical and, consequently the polarizability tensors are expected to be anisotropic. 
We also calculate the anisotropy in polarizability given by 



lAal 



3Tra 2 - (T 



rot 



1/2 



(general axes) 



(7) 



where a is the second-rank polarizability tensor. 

The two methods (MP2, and TDDFT), which we employ to calculate the static polariz- 
abilities of clusters (for n > 4) take into account the electron correlations in different ways 
and thereby enabling us to check the consistencies of our results, as no other theoretical 
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results are available in the literature for comparison. In Table V we present results for the 
average static polarizability a and anisotropy in the polarizability Aa obtained with MP2 
and TDDFT (with SAOP and LDA XC potentials) methods. We have performed MP2 
calculations of a and Aa for clusters only up to K 14 and K 10 , respectively. We note from 
Table V that the results for both polarizabilities and their anisotropies obtained by MP2 
and TDDFT-SAOP are quite close for all the clusters while corresponding LDA values are 
systematically underestimated. SAOP results for all the isomers for each cluster are actually 
slightly higher than the corresponding MP2 data except for K 2 . In order to study the evo- 
lution of polarizability a and anisotropy in polarizability Aa with the size, we plot them in 
Figs. 4a and 4b, respectively for the minimum energy isomers of each cluster as a function of 
number of atoms. In this figure we also display experimental data for three magic clusters, 
namely, K2, Kg, and K20 Fig. 4 once again clearly shows that MP2 results for the 
polarizabilites and anisotropies in polarizabilities are very close to the corresponding SAOP 
results and LDA values are systematically lower than both of them. Both SAOP and MP2 
results for the polarizabilities for magic number clusters K2, Kg and K20 are well within the 
experimental error bars (of the order of ± 5 — 7%). The maximum difference between the 
experimental and SAOP as well as MP2 results is observed for K 2 o cluster. Note also (from 
Table V) that both SAOP and MP2 results for the polarizabilities of clusters Kg and K20 
possessing T^ symmetry ( K 8 _! and K 2 o_4 ) are closer to the corresponding experimental data 
(for Kg, a = 1653 ± 83 a.u. and for K20 6t = 3834 ± 300 a.u. from experiments) than their 
respective minimum energy structures (see Table V). From these results it may be inferred 
that for these two clusters the geometries detected in the experiments performed at finite 
temperatures may be different from what have been obtained in this paper by DFT based 
calculations. Overall, we conclude from the results of Table V and Fig. 3 that TDDFT 
based calculations with SAOP yield results for polarizability which are quite accurate and 
compare well with the correlated wavefunction based MP2 results as well as experimental 
data. 

It has already been pointed out that MP2 calculations are computationally expensive 
and thus it becomes increasingly difficult to apply this method for very large clusters (more 
than 10 atoms). However, a good match between SAOP and MP2 results for clusters up to 
K 14 encourages us to explore how the two results scale with respect to each other. To this 
end we plot SAOP and MP2 results along y- and x-axes, respectively, in Fig. 5 and fit the 
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data points with a straight line by least square fitting. It can be clearly seen from Fig. 5 
that a very good fitting is obtained with the correlation coefficient value 0.9997 signifying a 
linear relationship between SAOP and MP2 results. This linear relationship between SAOP 
and MP2 will enable us to predict MP2 results for the polarizability of larger clusters. 

According to the iellium model, clusters with closed shells of delocalized electrons have 
spherica, shape B, i It is we» known that polarity of such a sphere is property 
to the volume of the sphere. Since the geometries of clusters considered in this paper are, 
in general, not spherical in nature and thus such linear dependence of polarizabiliy with 
the volume of cluster is not very obvious. However, the studies on the relationship between 
the static polarizability and the volume of carbon and sodium clusters have already been 



reported in the literature 



m 



15 



. Here we extend this study for potassium clusters 
and we go up to cluster containing 20 atoms. For this purpose we use SAOP results for the 
polarizabilities of lowest energy isomers for each cluster and obtain its volume by using the 
prescription of Tomasi and Persico 59J . The plot of the polarizability as a function of the 
volume of the clusters is shown in Fig. 6. It can be clearly seen from this figure that a good 
fitting is obtained with the correlation coefficient value of 0.996. This result then clearly 
suggests that a good correlation exists between the polarizability and the cluster volume 
even for nonspherical potassium clusters. This linear correlation between the polarizability 
and volume is an important result as it enables us to construct a size-to-property relationship 
for polarizability. Using this relationship polarizabilities of larger clusters can be calculated 
as for these clusters performing ab initio calculations are computationally expensive if not 
impossible. 

We wish to close this paper with a comparison of these SJBM results with the ab initio 
results presented in this paper. The jellium based results within DFT for the polarizability 



of K 8 , and K 20 are availab.e in the .iterate M and we — them with corresponding 
ab initio TDDFT results. The LDA calculations in Ref. [19j were performed by employing 
Dirac form for the exchange and Wigner form for the correlation energy functionals [60] . On 
the other hand, in this paper we employ VWN parametrization of the LDA XC functional 
which uses same Dirac exchange energy functional, but the parametrization for correlation 
part is different from the Wigner functional. We expect that the deviation in the results due 
to application of different correlation energy functionals will be significantly smaller than 
the difference in the two results, arising due to the consideration of structures of the clusters 
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in ab initio calculations. The results for the polarizabilities of Kg and K20 clusters calculated 
with LDA XC functional in jellium model are found to be ax 8 = 1212 a.u. and ax 20 = 2939 
a.u.. In comparison to this our ab initio results (for the minimum energy geometries) with 
LDA XC potential are a# 8 = 1367 a.u. and a^ 20 = 3220 a.u.. Besides LDA, Rubio et 
al. [l9j also employed a potential with correct — 1/r asymptotic decay as introduced by 
Przybylski and Borstel (PB) 6l| within weighed density approximation (WDA) to calculate 
the polarizability. We compare the results obtained with PB potential within jellium model 
( ®k$ = 1542 a.u. and «jf 20 = 3489 a.u.) with SAOP results (a^ g = 1531 a.u. and 
&k 2 q — 3582 a.u.) as SAOP too possesses correct asymptotic behaviour. The jellium based 
LDA results are around 10% lower than the corresponding TDDFT-based ab initio values. 
On the other hand, differences between ab initio SAOP and PB within jellium model results 
are further smaller. We note here that a similar observation for the polarizabilities of sodium 
clusters was made in Ref. llj. From the closeness of the results obtained by employing 



ab initio and SJBM, we conclude that for the alkali-metal atom clusters detailed ionic core 
structures may not have much influence on the values of the cluster polarizabilties. 



IV. CONCLUSION 



This paper is devoted to the calculations of the optimized geometries and various other 
properties of potassium clusters containing even number of atoms ranging from 2 to 20. In 
order to determine optimized geometry we have used DFT with TZ2P basis set and GGA 
XC potential. For each cluster beyond K 4 more than one geometry has been considered 
to explore the possibility of the existence of various structural isomers. To accomplish this 
we have taken various optimized geometries of sodium clusters available in the literature as 
our starting geometries. For small clusters containing 2, 4, and 6 atoms the results of our 
calculations for BE per atom, bond length, and average interatomic distance match quite 
well with the other published data obtained with correlated wavefunction based methods. 

We have also studied from these DFT based calculations evolution of IP with the size 
of clusters. The experimental data for the size dependence of IP of potassium clusters is 
available in the literature and our ab initio DFT based results match quite well with them. 

In this paper we have carried out calculations of the static polarizabiliteis of the potassium 
clusters with MP2 and TDDFT approaches as well as CCSD(T) (for dimer and tetramer) 
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taking electron correlations in different ways. A model XC potential (SAOP) possessing 
correct behaviours both in the asymptotic and inner regions of the molecule and also less 
accurate LDA XC potential have been used to calculate polarizabilities within TDDFT. For 
all the calculations sufficiently large basis sets have been employed. For dimer and tetramer 
the results for the polarizabilities obtained by different methods, employed in this paper, 
agree witheach other. Similarly for clusters beyond K 4 and up to K u MP2 and SAOP results 
for the polarizabilities are quite close to each other. We find a very good linear correlation 
between MP2 and SAOP results. On the other hand, TDDFT based calculations with LDA 
XC potentials are systematically lower than those of MP2 and SAOP. Moreover, we also 
find that both SAOP and MP2 results for the static polarizabilites of 2-, 8-, and 20-atom 
potassium clusters agree quite well with the experimental results. In general it is observed 
that both for sodium and potassium clusters the SAOP data for the polarizability is higher 
than the corresponding MP2 results. In this paper we have also investigated the volume- 
to-polarizability scaling for the potassium clusters. Our study has found a very good linear 
correlation between the volume and polarizability of the clusters. This scaling law can be 
exploited to determine the polarizabilities of larger clusters. Finally we have also compared 
the SJBM based results for the polarizabilities of K 8 and K 2 o with our corresponding ah 
initio values obtained by employing TDDFT. These comparison clearly reveals that jellium 
model based results for the polarizabilities are quite accurate for magic number clusters and 
it is expected that this model is increasingly more suitable for such larger clusters. 
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Figure captions 



Fig.l Optimized ground state geometries of K 2 , K 4 , and K 6 clusters. Letters denote the 
dimensions tabulated in Table I. 

Fig. 2 Optimized geometries of potassium clusters K n with n = 8 — 20. Binding energy 
per atom and average interatomic distance for these clusters are tabulated in Table II. 

Fig. 3 Plot of ionization potential (in eV) of potassium clusters containing even number 
of atoms ranging from 2 to 20 atoms as a function of number of atoms. The results obtained 
with BP86 XC potential ( solid squares) are compared with the experimental results (solid 
circles) and (solid triangles) of Refs. [18| and [48j respectively. The continuous dotted line 
shows the results obtained via Eq. [HI 

Fig. 4 Plot of average static polarizability (a) a and (b) anisotropy in polarizability Act 
for minimum energy isomers of potassium clusters as a function of number of particles. The 
experimental results for the polarizabilities of magic clusters K 2 , K 8 , and K 20 [1] are also 
shown in this figure. The lines joining the points are guide to the eye. 

Fig. 5 Plot of SAOP results for average static polarizability a against corresponding MP2 
values. All the results are in atomic units and straight line is least square fitted line. 

Fig. 6 Plot of average static polarizability a obtained with SAOP as a function of the 
cluster volume. All the results are in atomic units and straight line is least square fitted 
line. 
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TABLE I: Comparison of binding energy per atom ( in eV), bond length, and average interatomic 
distance (R) (in angstrom) for K2, K4, and Kq clusters. 



K 


Reference 


-BE/N 


Bond length 


(R) 


K 2 . 


Present 


0.25 


3.94 


3.94 




Ref. [29] 


0.19 


4.22 


4.22 




Ref. [27] 


0.21 


4.21 


4.21 




Ref. [28] 


0.32 


3.84 


3.84 




Ref. [30] 


0.26 


4.05 


4.05 




Ref. [31] 


0.27 


3.92 


3.92 




Ref. [45](Expt.) 


0.25 


3.90 


3.90 


K4.0 


Present 


0.31 


a= 4.44, b= 3.98 


4.34 




Ref. [29] 


0.28 


a = b = 4.78 


4.78 




Ref. [27] 


0.19 


a = b = 4.90 


4.90 




Ref. [28] 


0.37 


a = b = AAA 


4.44 




Ref. [30] 


0.34 


a= 4.42, b= 3.92 


4.44 


K 6 .o 


Present 


0.39 


a= 4.47, b= 4.31 


4.39 




Ref. [29] 


0.37 


a = b = 4.65 


4.65 




Ref. [28] 


0.44 


a = b = 4.50 


4.50 


K 6 .i 


Present 


0.39 


a= 4.29, b= 4.60 


4.39 




Ref. [29] 


0.36 


a = b = 4.65 


4.65 




Ref. [28] 


0.42 


a = b = 4.28 


4.28 
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TABLE II: Binding energy per atom ( in eV), and average interatomic distance (R) (in angstrom) 
of potassium clusters containing 8 to 20 atoms. 



K 


-BE/N 


(R) 


Kg.o 


0.458 


4.44 


K 8 .i 


0.448 


4.45 


Kio.o 


0.457 


4.48 


K10.1 


0.454 


4.50 


KlO.2 


0.453 


4.54 


Kio.3 


0.413 


4.50 


K12.0 


0.476 


4.54 


K12.1 


0.473 


4.50 


K12.2 


0.471 


4.51 


K14.0 


0.493 


4.57 


K14.1 


0.487 


4.50 


Kl4_2 


0.486 


4.48 


K16.0 


0.500 


4.57 


K16.1 


0.499 


4.57 


K16.2 


0.498 


4.53 


K18.0 


0.522 


4.52 


K18.1 


0.518 


4.60 


K18.2 


0.515 


4.62 


K20.0 


0.533 


4.53 


K20.1 


0.532 


4.57 


K20.2 


0.530 


4.59 


K20.3 


0.529 


4.52 


K20.4 


0.514 


4.56 
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TABLE III: Comparison of the theoretical and experimental IP (in eV) for K 2 , K4, and Kq clusters. 



Reference 


K 2 


K 4 


K 6 


Present (BP86) 


4.28 


3.61 


3.73 


Ref. [27] 


3.79 


3.23 




Ref. [28] 


4.12 


3.64 


3.79 


Ref. [30J 


4.32 


3.43 




Ref. [18] (Expt.) 


4.05 


3.52 


3.35 


Ref. [48] (Expt.) 


4.05± 0.05 


3.6± 0.1 


3.44± 0.1 


Ref. [46] (Expt.) 


4.05± 0.05 


3.6± 0.1 




Ref. [47] (Expt.) 




3.36 


3.25 



TABLE IV: Comparison of the average static polarizabilty a (in a.u.) for K2 and K4, clusters. 



Method (Reference) 


K 2 


K 4 


SAOP (Present) 


471.2 


994.1 


MP2 (Present) 


483.2 


977.8 


CCSD(T) (Present) 


510.2 


991.0 


CCSD(T) (Ref. [3JJ) 


486.4 




CCSD(T) (Ref. [56]) 


502.1 




Expt. (Ref. [18]) 


486.5 




Expt. (Ref. [57]) 


500±40 
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TABLE V: Average static polarizability a and anisotropy in polarizability Aa of various isomers 



of potassium clusters in atomic units. 





MP2 


SAOP 


LDA 


K n 




Arc 


LI 


Arv 


Lt 


Arv 


K 2 


483.2 


369.1 


471.02 


336.9 


437.4 


289.7 


K 4 


Q77 8 


931 1 


QQ4 1 


937 5 


Q07 3 


876 5 


K 6 .o 


1 ^08 ^ 


7^4 OQ 


1 ^91 1 


783 99 
( 00. zz 


1 9D3 1 


799 ^ 
( zz.o 


K 6 .i 


i/|iq 
141o.o 


y 1 1 .0 


1 /i /in n 

144U.U 


1(111 Q 


looU.U 


a/17 1 
y4 / .1 


K« n 


14Q9 9 


^01 n 


1 ^1 9 


390 94 


1 3fi7 5 


9fi7 Q 


K 8 .i 


±UZU.y 


9D fi 


1 fi4Q 
±U4:y 


91 8^i 
z 1.00 


1 48Q 4 


1 7 4^ 


Kio o 


ZUUO. 1 


1 OH/1 Q 

1ZU4.0 


ZUDZ.D 


1 1 7Q O 


io4y.4 


11Z0. / 


Km i 


9D1 ^ 4 


1 987 9 
1ZO ( .z 


90^9 8 
zuoz.o 


1 988 1 

1ZOO. 1 


1 84Q 8 


1 991 1 

1ZZ 1 . 1 


Kiq_2 


90^4 n 




908^ fi 


1 ifii n 


1 878 fi 
10 ( o.u 


111^1 


Kiq_3 


9^8 a 


fi8 e; 

Oo.O 


9/190 
Z4ZU.0 


on 7n 
yu. ( u 


9990 
ZZoo.o 


Q7 7 

y / . / 




9^8^ 4 




949Q 4 


1 4 


91 Q4 fi 


1 ^1 ^ 7 
1010. ( 


K12.1 


9448 9 
z^^o.z 




948^ 1 


1 4fi1 8 


99^Q fi 
zzoy.u 


loyu.o 


Ki2_2 


9^80 7 
ZooU. 1 




9/1 1 7 Q 

Z4i ( .y 


1 /Ifi7 8 
140 ( .0 


91 8fi 9 
ZloD.Z 


1 ^7^ 
10 ( 0.0 


*-±4t_U 


9781 fi 

Z ( Ol.U 




9789 
z ( oz.o 


1 ^1 fi 3 


zuou.^ 


1 4^4 


K14.1 


97fifi fi 




97fi9 ^ 


1 ^fi Q 


9^1 7 Q 


1448 Q 


Ki4_2 


Z001.0 




97QQ Q 
Z 1 tJO.y 


1 9^ fi 
1Z00.D 


9/179 8 
Z4 ( Z.o 


1 9n^ q 
izuo.y 


Ja.1 n 






988^ 9 


89^ 4 


9fi1 Q fi 


778 fi 
( ( o.u 


Kie.1 






9878 Q 
zo ( o.y 


81 fi Q 


9fi14 ^ 


7fiQ 4 


Kig_2 






9QQD Q 


Q1 9 4 


9799 3 
Z ( zz.o 


878 4 

0(0.4: 








^n/ia 7 
oU4o. 1 


/IS"? 1 
4oo. 1 


9787 7 
Z ( ( . ( 


/i^fi n 

400. u 


Kis.1 






391 fi 3 
OZIO.O 


R93 7 
OZo. ( 


9Q1 9 ^ 
zy iz.o 


OOo.U 


Ki8 2 






3197.7 


801.9 


2894.0 


744.4 


K 2 o.o 






3582.2 


344.40 


3220.4 


184.4 


K20.1 






3604.3 


140.1 


3289.1 


120.8 


K20.2 






3526.4 


539.3 


3203.0 


475.2 


K20.3 






3518.6 


204.4 


3220.4 


184.8 


K20.4 






4066.1 


6.35 


3769.8 


6.58 
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